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Abstract 

We study the interaction of light with two Bose condensates as an open quan- 
tum system. The two overlapping condensates occupy two different Zeeman 
sublevels and two driving light beams induce a coherent quantum tunneling 
between the condensates. We derive the master equation for the system. It is 
shown that stochastic simulations of the measurements of spontaneously scat- 
tered photons establish the relative phase between two Bose condensates, even 
though the condensates are initially in pure number states. These measure- 
ments are non-destructive for the condensates, because only light is scattered, 
but no atoms are removed from the system. Due to the macroscopic quantum 
interference the detection rate of photons depends substantially on the rela- 
tive phase between the condensates. This may provide a way to distinguish, 

whether the condensates are initially in number states or in coherent states. 
03.75.Fi,42.50.Vk,05.30.Jp,03.65.Bz 
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I. INTRODUCTION 



In the Bose-Einstein condensation (BEC) phase transition the Bose gas is expected to 
acquire nontrivial phase properties. In the conventional reasoning the condensate is given 
a macroscopic wave function which acts as a complex order parameter with an arbitrary 
but fixed phase jT],3. The selection of the phase implicates the spontaneous breakdown of 
the U(l) gauge symmetry, breaking the degeneracy of the ground state for a system. The 
condensate contains a macroscopic number of particles and the strict conservation of the 
number of atoms is abandoned. The removal or addition of one particle in the condensate is 
not assumed to affect essentially the system, resulting in non-vanishing expectation values 
for particle annihilation and creation operators in the ground state. These statements are 
rigorously valid only in the thermodynamical limit, where the number of particles N — > oo. 
However, real atom traps have a finite number of atoms, much less than the classical limit. 

Javanainen and Yoo |5) have shown that interfering two Bose condensates on an atom 
detector builds up a quantum coherence between the two condensates resulting in an in- 
terference pattern, even though the condensates are taken to be in number states with no 
phases and no violation of the particle number conservation. This was done by simulating a 
quantum measurement process for 1000 atoms. The resulting interference is a consequence 
of the correlations between atomic positions and the particular measurement process. It is 
necessary that the detected atoms from the two condensates are indistinguishable, i.e. one 
does not know which condensate the detected atom came from. This uncertainty of the rel- 
ative atom numbers between the condensates breaks the global gauge invariance associated 
with the atom number conservation. Recently, there has been number of related works on 
simulating numerically quantum measurements of the relative phase between two Bose con- 
densates PHlOU (see also [pl[] ) and photon outputs from single mode cavities [12|]. Perhaps, 
the most relevant to this paper is the consideration of the Josephson coupling between the 
two condensates ||. Despite the high scientific activity all the condensate phase simula- 
tions have been based on straightforward and rather idealized models of atom detection. In 
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particular, the measurements have been destructive, i.e. during the measurement process 
atoms have been removed from the condensates and the total number of atoms in the two 
condensates has not been constant. 

Several detection schemes of spontaneously broken gauge symmetry in BEC of atomic 
gases have been proposed. Javanainen |Tj| has shown that the amplification of phase coher- 
ent laser beams driving Raman transitions between two Bose condensates is a signature of 
the broken symmetry. Imamoglu and Kennedy |14| found that coherent spontaneous Raman 
scattering between two independent, spatially separated condensates may be eliminated by 
adjusting the phase difference of the driving laser beams. In our previous work it was shown 
that the relative peak heights in the spectrum of scattered light could possibly be used to 
determine the relative condensate phase JT5|, and that the spontaneous decay rate of inde- 



pendently produced excited atoms depends strongly on the spontaneously broken symmetry 
|16| . In all these papers the condensates were assumed to have macroscopic wave functions 



with well-defined phases in accordance with the spontaneous breakdown of gauge symmetry. 

In this paper we consider a gedanken experiment closely related to the previously pro- 
posed detection schemes of spontaneously broken gauge symmetry in BEC JT3|-|IT)|| . We 
evaluate the conditional probability for detecting coherent spontaneous scattering of pho- 
tons from two Bose condensates driven by external light beams. Spontaneous scattering 
specifies that the emission is not stimulated by light, although it is stimulated by atoms 
||16|| . The conditional probability depends on the macroscopic quantum coherence and it 
could possibly be used as a method to detect the phase difference of the two condensates. 
However, in our calculations no spontaneously broken gauge symmetry is initially assumed, 
but the relative phase of the two Bose condensates is established via direct measurements 
of spontaneously scattered photons. This measurement technique is non- destructive for the 
condensates; the light scatters atoms between the condensates, but the total number of atoms 
in the two condensates is conserved. Even though we start initially from number states with- 
out any phase information and without the violation of the particle number conservation, 
the condensates behave as if they had a phase. A non-destructive optical detection of a 



Bose condensate has been reported recently JT^|. Coherent forward scattering was measured 
by blocking the transmitted probe beam by a thin wire. Our scheme also relates closely to 
the experimental production of two overlapping condensates in the same trap using nearly 
lossless sympathetic cooling of one state via thermal contact with the other evaporatively 
cooled state |T8|] . The two species were the \F — l,m — —1) and \F — 2, m = 2) hyperfme 



spin states of Rb. Similarly to Refs. ||13l , |r5f we assume that the atoms are confined in the 



same trap and that they occupy two different Zeeman sublevels which are optically coupled 
through a common excited state by two low intensity off-resonant light beams. 

In the recent measurements Andrews et al. [|19[] have found evidence of macroscopic 
quantum coherence in BEC by interfering two condensates. The two condensates had the 
same internal atomic states, but different external wave functions. The spatially separated 
condensates were created by focusing a blue detuned laser into the center of the magnetic 
trap, generating a repulsive optical force. The interference fringes were measured optically 
by short resonant laser pulses (absorption imaging), after the condensates had expanded 
ballistically about 40 ms. If the condensates were initially in number states, the observed 
interference was a consequence of the uncertainty in which condensate the atoms absorbed 
photons. Similarily, in our scheme the macroscopic quantum interference results from the 
uncertainty in the initial state of atoms, when they are excited by the driving light beams. 

We begin in Sec. [J by introducing the Hamiltonian for the system. In the limit of large 
detunings of the driving light fields from the atomic resonances the excited state operators 
may be eliminated adiabatically. We obtain an effective two-state Hamiltonian between 
the two Bose condensates. The dynamics of the two condensates and the driving light is 



considered as an open quantum system in Sec. |T|. The many-particle master equation is 
derived by eliminating the vacuum modes of the electromagnetic fields. In Sec. [TV] we discuss 
the light-induced coherent quantum tunneling between the Bose condensates. In Sec. |V| the 
evolution of the master equation in terms of stochastic trajectories of state vectors is studied. 
We show that the conditional probability for the detection of the (n + l) th photon during 
the time interval [t,t + St], given that the n th detection occured at the time t n depends 



strongly on the relative phase between the condensates. The results of simulations are 
presented in Sec. [VT]. In the simulations the relative phase between the two Bose condensates 
is established by the measurements of spontaneously scattered photons, even though the 
condensates are initially in pure number states. We also propose a test which could possibly 
distinguish an initial number state from an initial coherent state, i.e. to test the validity of 
the conventional symmetry breaking arguments in BEC of dilute atomic gases. Finally, a 



few concluding remarks are made in Sec. |VI] . 



II. EFFECTIVE TWO-LEVEL HAMILTON! AN 

We consider two spatially overlapping Bose condensates, whose ground states differ in 
their internal quantum numbers. We assume the condensates in two different Zeeman sub- 
levels \b) = \g, m) and |c) = \g,m— 2) [p~3| , p~5j] , see Fig. [I]. The state c is optically coupled to 



the electronically excited state |e) = \e, m — 1) by the driving field E^ 2 having a polarization 
cr + and a dominant frequency ^2- Similarly, the state b is coupled to e by the driving field E^i 
with a polarization er_ and a dominant frequency Q\. We assume that the phase-coherent 
light fields E dl and E d2 propagate in the positive z direction and are detuned far from the 
resonances of the corresponding atomic transitions. The light fields are also assumed to 
be in coherent states. We only consider the coherent spontaneous scattering between the 
condensates, which is stimulated by a large number of atoms in the condensates. The decay 
into the non-condensate center-of-mass (cm.) states is also stimulated by the Bose-Einstein 
statistics. However, at very low temperatures this stimulation is much weaker, because most 
of the particles are in the condensates. The Bose stimulated spontaneous emission into the 
non-condensate modes also scatters photons into the solid angle of Air, whereas in the decay 
into the condensates the photons are emitted into a narrow cone in the forward direction. 
In addition to the Bose stimulation of the spontaneous emission there is the unstimulated 
free space decay 7, which is always present. This also scatters photons into the solid angle 
of An, and with a sufficiently large number of atoms in the two condensates the free space 
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decay may be ignored. The Hamiltonian for the system reads f[5,|20,21 



H = htUcb c'c + H u eb + £k) e k e k + ^ a l a g 



- £ ( / d 3 r d be ■ E 1 (r)0*(r)0 ek (r) 6^ + H.c 
" E (/ ^ d ce ■ E 2 (r)0:(r)0 ek (r) C te k + H.c.) , (1) 

where the creation operators for the atomic states e, b, and c are given by e\ tf, and c 1 , 
respectively. The condensate wave functions, which would typically be solutions of the 
Gross-Pitaevskii equation 0, are <p b and C . The excited state wave function for the cm. 
state k is ek with the cm. energy e k . The photon annihilation operator for the mode q is 
a q . The internal energies are described by the frequencies of the optical transitions e — > b 
and e — > c, which are u eb and c<j ec (c<j c fc = uo eb — u ec ) . The last two terms in Eq. (ffl) are for the 
atom-light dipole interaction. The dipole matrix element for the atomic transition e — > b is 
given by d be . 

In the limit of large detunings, Ai = Qi — u) eb and A 2 = fi 2 — ^ ec , the excited state 
operators e k in Eq. ([I]) may be eliminated adiabatically, and the cm. energies of the excited 
state may be ignored [OT. We insert the steady-state solutions of the Heisenberg equations 
of motion for the slowly varying operators e k = e* ni 'e k into the Hamiltonian (|I|) as in 
Refs. [[y|[l5!. With the help of the completeness of the cm. states ek , Eq. fll]) reduces to 



an effective two- state Hamiltonian 

H = hu cb + huJ q a\a q - — ^b J d 3 r d be ■ Ei(r) d eb ■ Ei(r)#J(r)0 6 (r) 
+ C t c | tfV d ce • E 2 (r) d ec • E 2 (r)#(r)0 c (r) 

+ (tfc J d 3 rd be ■ Ei(r) d ec • E 2 (r)0*(r)0 c (r) + H.c.)} . (2) 

The electric fields may be solved from Eq. (0) by integrating the Heisenberg equations of 
motions ||20|| . It is advantageous first to postpone the rotating wave approximation (RWA) 
and then to define the positive frequency parts of the light fields by the dominant time 
dependence of the operators. The positive frequency component of the total electric field 
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= + is expressed in terms of the classical driving field E^, that would prevail 
in the absence of matter, and the scattered field E^. The radiation reaction effects |2T| 



may be ignored in the limit of large detuning, when the multiple scattering is negligible. 
The driving field inside the sample in the length gauge should be understood as the driving 
electric displacement divided by the permittivity of the vacuum |2T|J2^] . Assuming the light 
is measured outside the sample the scattered fields are given by 

E+(r) = ^i- J dV K(d 6e ; r - r>*(r') {0 6 (r')d e6 • E+(r') tfb + c (r')d ec ■ E+ (r') tf~c} , 

(3a) 



X — J rfVK(d ce ;r-rO0:(rO{0 c (rOd ec -E+(rO£ t £ + fe (rOd e6 -E+(r')5 t 6} , 

(3b) 



where we have defined the slowly varying operators Ef = e^'Ef, E^ = e^^E^, and 
c = e i (^i-^2)* c h ave made the first Born approximation and replaced the electric 

fields under the integrals in Eq. (||]) by the corresponding driving fields. The kernel K(T>; r) 
coincides with the classical expression |2]| of the positive-frequency component of the electric 
field from a monochromatic dipole with the complex amplitude T>, given that the dipole 
resides at the origin and the field is observed at r ^ 0. The explicit expression is 

K(X>; r) = {k 2 (hxV) x n— + [3n(n • T>) — T>](— — ^)e* r }, (4) 

47re r r 6 r z 

where k = Q/c and n is a unit vector pointing from the source point toward the field point. 



III. MASTER EQUATION 

In this section we consider the dynamics of the two Bose condensates and the driving 
light fields as an open quantum system and eliminate the vacuum modes. We assume the 
driving electric fields to be plane waves 

E+(r) = i^e_ e — , E+ (r) = \S 2 e+e^ . (5) 



We insert = E£ + E£ and E^~ = Ej 2 + E+ 2 into the Hamiltonian (g), where E+ x 
and E^ 2 are now considered as vacuum fields, and keep only the terms of second order 
in 1/Ai, i.e. we neglect the terms, which contain a product of two scattered fields. The 
system Hamiltonian H$ for the condensates and the macroscopic light fields is coupled to 
the reservoir of the vacuum fields by the Hamiltonian Hsr. The total Hamiltonian is given 
by H = Hs + Hr + Hsr, where the different parts have the following expressions 

HSR = "7iA7 {/ ^ deb ' E ^ l(r) ( dbe ' E ^ l(r) ^* (r) ^ (r) 6+6 + dce • E ^( r )C( r )<M r ) Jb) + H. 

{ / d 3 r d ec ■ E+(r) (d ce ■ E d2 (r) #(r)0 e (r) c t c + d be • E^(r) fl(r)^(r) &t c ) + H .cj , 
ffs* = ^(ciJcfc — 62) c^c — M\b^b + (hn b^c + H.cj , Hr = frtUg a) q a q . 



Here we have used the following notation 

- S l^/^«( r )e--> c ( r ,, (7, 

where K12 = k,i — k 2 is the wavevector difference of the incoming light fields, and 5i and 52 
are the light-induced level shifts. The dipole matrix element d eb contains the reduced dipole 
matrix element and the corresponding nonvanishing Clebsch-Gordan coefficient. 

We assume a reservoir at T = and obtain the equation of motion for the reduced 
density matrix of the system in the interaction picture ps using the standard Born and 
Markov approximations 

Ps = J™ dTTrR{[H SR (t),[Hs R {t-T),p s {t)®p R ))} , (8) 

where the trace is calculated over the vacuum modes and the reduced density matrix of 
the reservoir p R is assumed to be in a statistical mixture of eigenstates of Hr ||24 |. The 
calculation of Eq. (d) reduces to the evaluation of the vacuum field correlations 

1 t°° 

P(r 12 ,n u d eb ) = -^ dTe^ T (d eb --E+(T 1 ,t)d he -Ej 1 (r 2 ,t-r)). (9) 
n Jo 

Lenz et al. |2H| have derived the master equation for light matter interactions in the RWA. 



Applying the RWA leads to spatially non-local Cauchy principal value terms, which also 
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come out in the integration of the scattered electric fields (f|), in addition to the dipole 
radiation kernel, if the straightforward RWA is performed |[20|| . To avoid obvious physical 
and mathematical problems of the RWA we include in the vacuum electromagnetic fields 
the negative mode frequencies (u q < 0) for the photon annihilation and creation operators. 
This is basically equivalent to ignoring some of the commutation relations between photon 
annihilation and creation operators. 

We assume that there is a cutoff in the wave numbers q of the photons; we multiply the 
density of the states of the electromagnetic fields by e~ q2a2 ^ A , with a > being a length 
scale. The cutoff removes all mathematical problems concerning, e.g., the exchange of the 
order of derivatives and integrals, which are abundant in the theory without the cutoff. At 
the end of the calculations we ultimately take the limit a — > 0. 

The evaluation of Eq. (|9]) is very similiar to the integration of the electric field in 
Ref. |20| , |21|j . The time evolution of the vacuum fields in the interaction picture is defined by 
Hr. First, we set r 12 = |r^ — r 2 | = 0. For a small but nonzero a, the result is 

P(r 12 = 0, n 1? d eb ) = -^L + z^|^3 ee leb + i 8u eb (10) 

The first term ^ eb in Eq. (|10"D is the familiar spontaneous linewidth of the atomic transition 
e — > b. The second term Su eb diverges as the photon momentum cutoff goes to infinity with 
a —>■ 0. This part, after a proper renormalization, contributes to the Lamb shift. From now 
on we assume that the Lamb shifts are already included in the transition frequencies, and 
ignore the 5u eb term in Eq. (|T0|). 

The expression (|) is not divergent for r i2 ^ 0. In that case we obtain 

ic e ifc i r 'i2 ^ 

P(n 2 , fix, d eb ) = — d eb ■ (d be x V) x V = -- d eb ■ K(d be ; r 12 ) , (11) 

Awe^n ri2 n 

The value of not having made the RWA in the evaluation of Eq. (§) emerges in the apparent 

physical results Eqs. ([TOD and (JTTJ) . 

By defining the operators 

(7 x (r) = (d eb ■ E+(r) ^(r)0 6 (r) b% + d ec • E+ (r)0*(r)0 c (r) fetg) , (12a) 
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C 2 (r) = (d ec • E+ (r) #(r)0 e (r) e*c + d eb • E+(r) #(r)&(r) g f &) , (12b) 

the equation of motion for the reduced density matrix is obtained from Eqs. (Q), (|K|, 
and (0) 

z ~ 

+ i | rf 3 n d 3 r 2 d eb ■ K(d eb ; ria) (C 1 t (r 1 )C 1 (r 2 )p 5 - C 1 (r 2 )p 5 ^i t (ri) 
_ * | d 3 n rf 3 r2 ^ . K * (cU ria) ^^(rOCiCra) - C 1 (r 2 )p 5 C 1 t (r 1 
+ LJ d 3 ri d 3 r 2 d ec • K(d ec ; r 12 ) (C 2 t (r 1 )C 2 (r 2 )p 5 - C 2 (r 2 )p 5 C 2 t (r 1 ) 
_ * y rf 3 ri rf s r2 ^ . K * (dec . ri2) (p 5 C 2 t( ri )C 2 ( r2 ) _ C 2 (r 2 )p 5 C 2 t (r 1 )) (13) 

where the terms containing K(d; ri 2 ) should be understood in such a way that the divergent 
in-phase part of the dipole field at r 12 = is ignored. In the limit r 12 — > the real part 
of iK(d e {>; ri 2 )/7l coincides with 7 e 6, the real part of Eq. ([10|). The Hamiltonian ifs is the 
system Hamiltonian Eq. (|6]) in terms of the slowly varying operators c 

# 5 = -h(5 cb + 5 2 ) c ] c - hSrfb + (hK 6 f c + H.c.) , (14) 

where we have defined the two-photon detuning 5 c b = fl\ — fl 2 — u cb . 

The non-local interaction terms in Eq. fll~3|) depending on K(d; r 12 ) represent the dipole- 



dipole interactions between atoms. According to Ref. | 20| , if the detunings Ai and A 2 
are at least comparable to the collective linewidth of the Bose condensate T, and if the 
characteristic size of the condensate I is substantially larger than the wavelength of light 
I 3> A, the dipole-dipole interactions and the multiple scattering of light may be ignored. 
The condensate collective linewidth has been estimated to be T = 3N / y/(2l 2 k 2 ), where 7 is 



the transition's free space natural linewidth and N is the total number of atoms [^,^]. We 
assume these conditions to be valid and neglect the dipole-dipole interactions between the 
atoms. Because in the limit of low light intensity the contact interactions are inconsequential 
in the interactions between dipole atoms P2L this is done by ignoring the terms depending 



on K(d; r 12 ) with r 12 7^ in Eq. (0). We obtain for the master equation 
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Ps 



~[Hs, Ps] ~ leb J d 3 r (C 1 t (r)C 1 (r)p 5 + p^r^r) - 2C 1 (r)p 5 C 1 t (r)) 
- lec J d?r (pl{v)C 2 {v) Ps + psCl(v)C 2 (v) - 2C 2 (r)p 5 C 2 t (r)) . (15) 



IV. JOSEPHSON EFFECT 

In the absence of the coupling to the reservoir the dynamics of the system determined 
by the system Hamiltonian H$ in Eq. (|6|) can be solved analytically. The solutions for the 
ground state annihilation operators Ck and 6k are 

c(t) = e *(*+*i)* |c(0) ^cos Q R t + ^- sin Sl R tJ - ^-6(0) sin fi fl t | , (16a) 

&(t) = e *(*+*iX J 6(0) (cos Q R t - ^- sin n R t] - ^c(O) sin 1 , (16b) 

where we have defined "one-half of the effective two-photon detuning" <5 = (5 c f> — 5i +^)/2. 
The oscillation frequency in Eq. (O) is given by Cl R = (5 2 + k 2 ) 1 ^ 2 . To simplify the algebra, 
we have assumed k to be real. 

Before the light is switched on, the atoms in the states b and c are assumed to be 
uncorrelated. The driving light fields induce a coupling between the two Bose condensates, 



which is analogous to the coherent tunneling of Cooper pairs in a Josephson junction [28 
According to the Josephson effect, the atom numbers of the condensates oscillate even if 
the number of atoms in each condensate is initially equal. This can be seen easily by using 
the basic spontaneous symmetry breaking arguments and assuming the condensates to be 
in coherent states. Then, the annihilation operators have nonvanishing expectation values 
at the time t = 0: 



(6(0)) = ^jN b e^\ (c(0)> = ^jN c e^ . (17) 

We have assumed that the condensates in the states b and c have the expectation values 
for the number operators Nf, and N c , respectively. Now, if Nj, — N c — N/2, where N is the 
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total number of atoms, the expectation value for the occupation number in the condensate 
c is obtained from Eq. fllED 

(&(t)c(t)) = — f 1 - — sin {ip) sin {2Q R t) - ^- cos (<p) sin 2 (Q B t)j , (18) 

where if = ip c — ip b . The oscillation in Eq. (|T^) is a consequence of the macroscopic quantum 
coherence. In each measurement process the relative phase between the two condensates is 
selected as a random number. The value of the phase difference determines the amplitude 
of the oscillations according to Eq. ([18|) . 



V. STOCHASTIC SCHRODINGER EQUATION 



In this section we study the evolution of the master equation ([15]) in terms of stochastic 
trajectories of state vectors p9|-|3T|. The dissipation of energy from the quantum system of 



macroscopic light fields and the two Bose condensates is described by the coupling to a zero 
temperature reservoir of vacuum modes, resulting in the spontaneous linewidth for atoms. 
The master equation (115]) is in the Lindblad form |32] and is equivalent to the Monte-Carlo 



evolution of wave functions (MCWF) . The MCWF procedure consists of the evolution of 
the system with a non-hermitian Hamiltonian H e g, and randomly decided quantum 'jumps', 
followed by wave function normalization. In our case the quantum jumps correspond to the 
detections of spontaneously emitted photons. The system evolution is thus conditioned on 
the outcome of a measurement. The non-hermitian Hamiltonian is obtained from Eq. ( |15| ) 

H cS = H s -ih( leb J d 3 rC 1 t (r)C 1 (r)+ 7ec | d 3 r Cj(r)C7 2 (r)) . (19) 

The non-unitary evolution corresponds to the modification of the state of the system associ- 
ated with a zero detection result of the spontaneously emitted photons. Because the output 
is being continuously monitored by the detector, we gain information about the system even 
if no photons have been emitted. 

The operators C\ and C2, from Eq. flT2"D, correspond to the excitations of atoms by the 
driving light beams followed by the emissions of the photons with the polarizations <r_ and 
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cr + , respectively. In the translationally invariant system and for a non-interacting gas the 
spontaneous emission, Eq. (Q), is directed exactly parallel to the positive z axis, because of 
the momentum conservation. In that direction the two polarizations are perfectly distin- 
guishable. In a finite size trap the uncertainty in momenta introduces a narrow scattering 
cone for spontaneously emitted photons. For simplicity, we assume that this cone is narrow 
enough, so that the two polarizations can still be approximated as perfectly distinguishable. 
We also assume that all the spontaneously emitted photons can be detected. Then, we have 
two detection channels corresponding to the two different polarizations. 

The state vector at the time t is denoted by ip S ys{t)- Its evolution is determined by 
the non-hermitian Hamiltonian H c q, defined in Eq. (p!9[). If the wave function ip sys (t) is 
normalized, the probability that a photon with the polarization <r_ is detected during the 
time interval [t, t + 5t] is 

P_ = 2 leb J d 3 r (^ y s(t) I C 1 t (r)C 1 (r) | ^ ys (i)) 8t . (20) 

Similarly, the probability for the detection of a photon with the polarization cr + is 

P + = 2 lec J d 3 r^ sys (t)\Cl(r)C 2 (r)\^ sys (t))5t. (21) 

The probability of no detections is 1 — P_ — P + . 

In the implementation of the simulation algorithm we first, at the time t , generate a 
quasi-random number e which is uniformly distributed between and 1 . We assume that the 
state vector f/Vs^o) a t the time to is normalized. Then, we evolve the state vector by H e g 
iteratively for finite time steps At ~ 5t. At each time step n we compare e to the reduced 
norm 

(^ sys (t + nAt) | ^ sys (t + nAt)) = (^ sys (t ) I e iH > M / h e- iH °« nAt / n | ^ sys (t )> > ( 22 ) 

until (ipsysito + t) I V'sys^o + 1)) < e, when the detection of a photon occurs. It is easy to 
show that, at first order in At, Eq. (^) is the joint probability of not detecting photons in 
any of the n time intervals. After the detection we generate a new quasi-random number 
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rj. We evaluate P_ from Eq. fl2"0|) and P + from Eq. ( pT|) at the time of the detection. If 
77 < P_/ (P_ + P+), we say the polarization of the detected photon is cr_. If the photon has 
been observed during the time step t — > t + At, we take the new wave function at t + At: 

I ij sys (t + At)} = ^ h I d 3 r d(r) I Vsys(t)} , (23) 



which is then normalized. Otherwise, rj > P_j (P_ + P+) and the polarization of the detected 
photon is cr + . In that case the new wave function before the normalization reads 

I Vs ys (t + At)) = f d 3 rC 2 {r) | ^ sys (t)) . (24) 



After each detection the process starts again from the beginning. 

The conditional probability for the detection of the (n + l) th photon during the time 
interval [t + t n ,t + t n + St], given that the n th detection occured at the time t n , follows from 
Eqs. (H), ©, and © 



Pit + t n , t n ) = P_(t + £ n ; t n ) + P + (t + t n , t n ) , (25) 
where the expression for P_(i + £ n ; in) is given explicitly by 
P_(t + t n ;t n ) 

= 2/y b f d 3 r (Ays(t n ) I exp {iHl s t/h} Cl(r)d(r) exp {-iP cff t/ft} | i/j ByB (t n )) gf 
(4> sys {t n )\ exp{iHl s t/h}exp{-iH cS t/h}\4} sys (t n )) 

~ 2 7efe f d 3 r (Vsys(tn) I e-^ st " /n Ct(r, t + t n )d(r, t + t n ) e l ^ n \ ^ sya (t n )) St . (26b) 



Here we have written the operator Ci(rt) = exp {iH$t /fi}C i(r) exp {— iHgt/h} in the 
Heisenberg picture with respect to the system Hamiltonian H$, defined in Eq. fll4|). Equa- 
tion ( 26b ) is valid, if the characteristic time evolution of the system Hamiltonian is much 



faster than the detection rate of the spontaneously emitted photons |J . The expression (^6 
is especially useful, if the measurements of the spontaneously scattered photons Eqs. (|23| 
and (01) do not significantly disturb the system, i.e. if 



e i Hs t n+1 /n 1 ^ sys(tn+1 )) ~ e iH 3 t n /n 1 ^ sys ( tn )) 
14 



Then, the dynamics of the conditional probability P_(*; t n ) for detecting a photon with the 
polarizability cr_ during the time interval [*, * + 5t] after detecting the n th photon at the 
time t n is determined by the system evolution of Cjf (rt)Ci(rt). This system evolution may 
be solved analytically from Eq. fll6]) . We assume that the expectation values for the number 
operators in both condensates are initially equal Nf, = N c = N/2, and that N 1. Then, 
the expectation value (C}(ri)(7i(ri))sB i n the spontaneous symmetry breaking state, where 
the condensate annihilation operators have nonvanishing expectation values Eq. fll7D , and 
the expectation value (C}(r*)Ci(r*))N in the pure number state are given by 

N 2 

(C 1 t (rt)C 1 (rt)) SB = — |4w(r) [1 + sin (2k*) sin (^)f 
lb l 

+v(r) [3 + cos (4k*) — cos (4k*) cos (2ip) + cos (2<p)] 

+8s(r) [1 + sin (2k*) sin (ip)\ cos (</?)} , (27a) 
N 2 

(C 1 t (r*)Ci(r*)) N ~ — {2w(r) [2 + sin 2 (2k*)] + v(r) [3 + cos (4k*)]} . (27b) 
Here we have used the following definitions 



w(r) = \<p(r)\^, v(r) = |0(r)| 4 ^, s(r) = ^|^|0(r)| . (28) 



To simplify the expressions, we have set in Eqs. fl27|) and (p8|) 5 = 0, 0;,(r) = c (r) = 0(r), 
and s(r) real. We can simplify Eq. ( p?D further by setting w(r) = v (r) = s(r) 

(C 1 t (r*)C 1 (r*)) SB = N 2 w(r) cos 2 f-) [1 + sin (2k*) sin (<p)\ , (29a) 



(C{(rt)Ci(rt)> N£ -— ii. (29b) 



The effects of macroscopic quantum coherence in Eqs. ( |27a|) and ( |29a| ) are clearly ob 



served, as they are in the Josephson coupling, in Eq. (|T8"D. If the relative phase between 
the two condensates <p = it, the expectation value in Eq. (|29a| ) completely vanishes for the 
perfect spatial overlap of the condensate wave functions. Then, the conditional probability 
P-{t;t n ) in Eq. Q26bp for detecting a photon with the polarizability cr_ during the time 
interval [*, * + St] after detecting the n th photon at the time * n also vanishes. If tp = 0, 
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the expectation value (|29a|) is time-independent and twice as much as the expectation value 
( [29 b| ) for the number state. For simplicity, we set 5 = in Eqs. (pTD, (f28|), and fl2£f). These 
expressions are valid provided that the frequency of the detections is much larger than 5. 
If the detuning from the two photon resonance becomes small, the effective linewidth 7 of 
the transition c — > b may have an effect. However, it may be shown to be proportional to 
A 1 or smaller. The same expressions for C\C 2 may be obtained from Eq. (|29|) by changing 
ip -(p. 

Because in the scattering processes of atoms between the condensates the initial states 
are indistinguishable, i.e. we do not know which condensate the atoms are scattered from, 
the transition amplitudes interfere. For a large number of atoms these interference effects 
may only be observed, if there is a macroscopic quantum coherence of atoms present in the 
initial states, as in Eqs. (|27aD and ( |29a|) . In Ref. |15| we investigated incoherent scattering, 
in which case there is an uncertainty in both the initial and the final states of atoms. The 
macroscopic quantum interference of the transition amplitudes affects both the scatterings 
of atoms into the condensates and out of the condensates. 

In this section we have shown that the macroscopic quantum coherence has clearly ob- 
servable effects on the conditional probability P(t; t n ) in Eq. for spontaneous emission 
of photons. What remains to be shown is that the measurement process of the sponta- 
neously scattered photons duly reaches a steady-state in which it does not significantly 
disturb the dynamics of the system, and that this steady-state coincides with the theory of 
the spontaneous breakdown of the gauge symmetry in BEC. 



VI. RESULTS OF SIMULATIONS 

We have simulated the light matter dynamics using the algorithms described in Sec. [V]. 
The two condensates are assumed to be initially in the number states with equal atom 
numbers, Nf, — N c — N/2. We define the visibility of the macroscopic quantum interference 
j3 and the angle <p as the modulus and the phase of the complex number after n detections 
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fie* = ^(^s ys (tn) I e-^' n &+(0)c(0) e i6 ^ h \ ^ sys (t n )) . (30) 

According to the spontaneous symmetry breaking arguments Eq. (|T7j), for a coherent state 
we have (3=1, and = ip c — ipb is the relative macroscopic phase between the two con- 
densates. However, for number states there is no phase information at all and (3 = before 
any detections are made. Starting from pure number states the simulations show that a 
macroscopic quantum coherence is established by measurements of spontaneously scattered 
photons. In each run of measurements the relative phase between the condensates is se- 
lected as a random number. Although the symmetry is broken in each individual run of 
measurements, it is regained when an ensemble of measurements is considered. 

With very small values of the effective two-photon detuning the system approaches a 
dark state and the time between the photon detections increases rapidly. In the dark state 
the driving light does not excite the atoms. The coherent quantum tunneling between the 
condensates Eq. flIE|) dies out and the conditional probability of detecting spontaneously 
emitted photons, Eqs. ( 2515 ) and ( |29a| ), goes to zero. This corresponds to a vanishing 



relative phase between the condensates. On the other hand, if the effective two-photon 
detuning is very large, the measurements drive the system into a state where most of the 
particles are in one of the two condensates. The maximum visibility of the interference 
between the condensates is reduced from one to /3 max = 2VNbN c /N, if the condensates 
have unequal atom numbers. According to Bose-Einstein statistics, the scattering to an 
already occupied state is enhanced. If most of the atoms are in the condensate b, the 
probability P of detecting a photon with the polarization o~_ Eq. (|20|) is much larger 
than P + , the probability of detecting a photon with the polarization cr + Eq. (PI|). For 



Nb ^> N c , we also have (C\Ci) ~ (b^bb^b). Thus, the measurements end up driving the 
system towards a number state reducing the macroscopic coherence. However, even with 
these two extreme values of the effective two-photon detuning the measurements build up 
a macroscopic quantum coherence with a randomly-defined relative phase between the two 
condensates, before the coherence starts decreasing or the original phase is lost. This is an 
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example of the similiar behaviour of an initial coherent state and an initial number state 
after a large number of measurements. The detections rapidly build up a coherence and a 
random phase value for the initial number state, even though these are not preserved for an 
arbitrary number of measurements. In the absence of collisions and with a carefully chosen 
two-photon detuning the visibility (3 ~ 1 may be maintained in the simulations practically 
for an arbitrary number of detections. 

We have run the simulations for the total number of 1000 atoms with the various values 
of two-photon detuning 5 C &. The light-induced level shifts 8\ and 62 from Eq. (0) are chosen 
equal to the coupling coefficient k, which is approximately 5 x 10 6 times the conditional 
probability of detecting the scattered photons from the pure number state iV 6 = N c in 
Eq. (|25|) . We have assumed that the external wave functions of the condensates are equal 
0&O) = c (r). 

In Fig. |2| we have shown the results of simulations of one run of 1000 detections. The 
two-photon detuning is given by 6 c b = 0.15k. The visibility f3 in Fig. 0(a) approaches 
to one very rapidly as the number of detections is increased. The small decrease in the 
visibility after 700 detections is a consequence of unequal atom numbers in the condensates, 
as the maximum visibility is reduced from one to /3 max . In Fig- @(b) /3sinyj also becomes 
well-defined as more detections are made. 

In Fig. |3] we have an example where the two-level detuning is larger, 5 c b = 2.0k. The 
visibility (3 in Fig. |](a) is close to one after 100 detections, but it starts then slowly decreasing 
as more atoms enter the Zeeman state b. The maximum visibility of the interference between 
the condensates is reduced, because of the unequal atom numbers. In Fig. ||(b) we have 
plotted the relative visibility defined by (3 r = /3//3 max . This remains close to one. The 
relative phase between the condensates in Fig. ^|(c) is reasonably steady. 

The results of simulations for the two-level detuning 5 c b = 0.05k are shown in Fig. |j. The 
system approaches a dark state due to spontaneous emission. The time elapsed from the 
beginning of the simulation increases in Fig. ^(a) very rapidly as a function of the number of 
measurements. According to Eqs. ([25|), (p7a|) , and (|29a|) the vanishing expectation value for 



the conditional probability of observing photons corresponds to a vanishing relative phase 
between the condensates. The value of the phase (p in Fig. |](b) rapidly approaches zero as the 
detection time starts increasing, although the visibility (3 remains close to one in Fig. f|(c). 
This is an example of the strong dependence of the state of the system on the measurement 
scheme. Although the measurements first define a random phase, the final value of the phase 
is not a random number. The detections of the spontaneously emitted photons drive the 
system into a state in which the relative phase between the two condensates vanishes. 

The conditional probability for the detection of the {n + l) th photon during the time 
interval [t, t + St], given that the n th detection occured at the time t n , Eq. (|25D , does not 



have a dependence on the relative phase of the two condensates, as in Eqs. ( |27a|) and 
( P9a|) , before the macroscopic quantum coherence is established by measurements. According 
to the conventional symmetry breaking arguments such a phase dependence should exist 
immediately from the first measurement. Given that the time between photon detections 
could be determined very accurately, a different behaviour would be observed, whether the 
condensates are initially in number states or coherent states. Similarly, if we start from pure 
number states, there is no coherent quantum tunneling between the condensates according 
to Eq. QTBP before p ^ 0. 

According to Eq. (|29| ) the conditional probability P(t; t n ) for detecting spontaneously 
emitted photons can completely vanish for the spontaneous symmetry breaking state, while 
it is always finite for a number state. If the relative phase between the condensates is equal 
to zero, no photons in the spontaneous symmetry breaking state are observed for perfect 
spatial overlap of the condensate wave functions. However, if the condensates were initially 
in number states, emitted photons would be detected according to Eq. ( |29b|) before the 
macroscopic coherence is established by measurements, even if this coherence suppressed 
the emission of photons. In our simulations for a thousand atoms the phase becomes well- 
defined after approximately 100-200 detections. With the total number of atoms N < 1000, 
the required number of measurements seems to be well above y/N. For condensates formed 



with from 10 6 to 10 s atoms, a significantly different behaviour is observed, whether the 
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condensates are initially in number states or in coherent states. This is an alternative 



scheme for testing Bose broken symmetry arguments to that suggested by Wong et al. [33 



They proposed that the validity of the Bose-broken symmetry could be tested from the 
collapses and revivals of the macroscopic quantum coherence [51, if atoms are not lost 
from the system during the measurement process. 

For our light scattering scheme to be a practical test of the validity of the spontaneous 
symmetry breaking arguments for dilute atomic gases, the spatial overlap of the condensate 
wavef unctions should be significant and the phase diffusion by collisions between different 
atoms should be slow compared to the spontaneous emission rate. According to Eq. ([29|) the 
probability of detecting spontaneously emitted photons as a function of time is stimulated 
by the large number of atoms in the condensates. Depending on the geometry of the system 
the spontaneous emission rate is roughly proportional to ^NgN ~ 7-/V 2 5i/Ai, for Si ~ 5%. 
Here N is the total number of atoms in the condensates, 7 is the free space spontaneous 
emission rate, N e is the number of electronically excited atoms, and Si and S2 are the light 
induced level-shifts from Eq. (0). The detuning of the incoming light Ai was chosen large. 
Wong et al. J7| have studied the effect of collisions on the relative phase between two Bose 
condensates. The two-body collision rate can be estimated by k ~ pira 2 v rms , where a is 
the scattering length, p is the density of atoms, and t> rms is the root-mean-square speed of 
atoms. At a temperature of 180 nK and with a density of 10 12 cm -3 this gives for 87 Rb 
atoms a collision rate of about one collision per second. In Ref. the simulations for 200 
atoms show a recognizable steady conditional phase for the collision rate k = 0.2r a , where 
T a is the atom detection rate. A strong overlap between the condensates could significantly 
reduce the decoherence effects of the collisions. 



VII. CONCLUSIONS 

We have studied carefully the interaction of light with two Bose condensates as an open 
quantum system including the effects of measurements. In the limit of large detuning of 
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the driving light beams from the atomic resonances, we have shown that the relative phase 
between the two Bose condensates may be established by the measurements of spontaneously 
scattered photons, even though the condensates are initially in pure number states. In 
the quantum trajectory simulations the strong effect of measurements on the state of the 
condensates is clearly observed. Even the modification of the effective two-photon detuning 
changes substantially the effects of the measurements. 

The measurement of spontaneously scattered photons is non-destructive for the conden- 
sates, because only light is scattered, but atoms are not removed from the two condensates. 
In particular, a non-destructive detection allows repetitions of independent runs of measure- 
ments for the same condensates at different times. Establishing the relative phase between 
the condensates by the simulations of measurements on scattered light has also other ad- 
vantages over the simulations of atom counting. In the case of light scattering we can use 



the well-known theories of photon detection [p5[ . The present authors are not aware of the 



existence of similiar theories for atom detection. The driving light also introduces naturally 
the high frequencies for the system dynamics, which are required for the Markov and Born 
approximations in the derivation of the stochastic Schrodinger equations. Thus, we have 
obtained the first evidence of the macroscopic coherence properties of Bose condensates, 
initially in number states, by simulations which may be justified by well-established phys- 
ical theories. Although the simulations of atom counting have previously shown coherence 
properties, it was not evident a priori that the photon detections also should, because the 
different measurements may affect the system in a very different way. 

We have shown that the conditional probability of detecting spontaneously scattered 
photons as a function of time depends strongly on the relative phase between the two Bose 
condensates. This may provide a method to detect the relative phase and to give an unambi- 
gious signature of the macroscopic quantum coherence in BEC. The significant dependence 
of the spontaneous emission rate on the condensate phase difference could possibly also be 
used as a test of the spontaneous symmetry breaking of the global gauge invariance in dilute 
atomic gases, i.e. as a way to determine the true quantum state of the Bose condensate. 
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FIGURES 

FIG. 1. The level scheme of the system. Two Bose condensates are in two different Zeeman 
sublevels \b) = \g, rri) and |c) = \g,m — 2). The states |c) and \b) are optically coupled to the 
center-of-mass manifold of the electronically excited state |e) = \e,m — 1) by the far-off-resonant 
driving fields E^ 2 and E^i having polarizations a + and <r_, respectively. 

FIG. 2. One run of simulations for 1000 atoms, (a) shows how the visibility of the macroscopic 
quantum interference (3 approaches one as a function of the number of measured photons. In (b) 
we have plotted (3 sin ip as a function of the number of detections. The relative phase between the 
two Bose condensates ip also becomes well-defined. Here we have set S\ = 82 = k and the two-level 
detuning 8 c b = 0.15k. 

FIG. 3. One run of simulations for 1000 atoms with S± = 82 = k and the two-level detuning 
8 c b = 2.0k. In (a) the visibility (3 is close to one after 100 detections of spontaneously scattered 
photons, but starts then decreasing, because of the unequal number of atoms in the two condensates. 
However, the relative visibility (3 r in (b) remains close to one. The relative phase between the two 
condensates is plotted in (c). 

FIG. 4. One run of simulations for 1000 atoms with 8± = 82 = k and the two-level detuning 
8 c b = 0.05k. The time elapsed from the beginning of the simulation as a function of the detected 
photons in (a) starts increasing as the system approaches a dark state. The value of the phase ip 
in (b) approaches rapidly zero as the detection time starts increasing, although the visibility (3 in 
(c) remains close to one. 
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Fig. 1 (Ruostekoski and Walls) 
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